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We construct a variational wave function to study whether a fully polarized Fermi sea of ultracold 
atoms is energetically stable against a single spin flip. Our variational wavefunction contains short- 
range correlations at least to the same level as Gutzwiller's projected wavefunction. For the Hubbard 
lattice model and the continuum model with pure repulsive interaction, we show a fully polarized 
Fermi sea is generally unstable even for infinite repulsive strength. By contrast for a resonance 
model, the ferromagnetic state is possible if the s-wave scattering length is positive and sufficiently 
large and the system is prepared to be orthogonal to molecular bound state. However, we can not 
rule out the possibility that more exotic correlation can destabilize the ferromagnetic state. 



Whether a fcrmion system with repulsive interaction 
will become ferromagnetic is a long-standing problem in 
condensed matter physics. Early in 1930s', Stoner used a 
simple mean-field theory to predict that ferromagnetism 
will always take place with sufficient large repulsive in- 
teraction However, this conclusion is later challenged 
by Gutzwiller who took the short-range correlation into 
account 0. So far, except a few specific cases d, [3], 
there is no conclusive result on itinerant ferromagnetism. 
Recently, MIT group reports an experiment on itiner- 
ant fermions in an ultracold Fermi gas with large posi- 
tive scattering length close to a Feshbach resonance Q, 
and they attribute their observations to Stoner ferromag- 
netism by comparing to theories However, these 
theories are basically mean-field theory or a second-order 
perturbation, which neither include the Gutzwiller type 
short-range correlation nor consider the unitary limited 
interaction nearby Feshbach resonance. Moreover many 
of the experimental signatures can be reprodiiced qual- 
itatively by a non-magnetic correlated state [8|. Thus, 
it calls for a serious study including the effects of both 
correlation and unitarity in this problem. 

In this Rapid Communication we address the question 
whether a fully magnetized state is stable against a single 
spin flip. We compare the energy of -I- 1 spin-up par- 
ticles with that of one spin-down particle and N spin-up 
particles. A fully magnetized ferromagnetic state is def- 
initely unstable if we can find a variational state of the 
latter whose energy is lower. Similar idea has been used 
previously in studying the stability of Nagaoka ferromag- 
netism in the Hubbard model [1, |9|-[T]| , and attractively 
interacting Fermi gases with large population imbalance 
[13, ■ In this work, we will explore different realiza- 
tions of "repulsive interactions" in ultracold Fermi gases: 

(I) Single-band Hubbard model in a two-dimensional 
(2D) square or three-dimensional (3D) cubic lattice. The 



Hamiltonian H = Ht + H;„t, Ht 



(ij) 



and 



Hint = 



h.c, where (ij) are nearest-neighbor sites 
UE^n^tn^i {U > 0) . 

(II) Continuum model with finite-range interaction 



potential in three dimension. H = X^ko- ^kc|(.^Ck(T + 
(1/^^) Eq,k,k' V{k - k')4+ktCq-kiCq-k';Cq+k't- Hcrc 
Ek = k'^/{2m), Q is system volume, and l^(k) is the 
Fourier transformation of real space interaction poten- 
tial V{r). Let To be the interaction range, and thus 
F(k) -> for |k| > fco = 1/ro. For the conve- 
nience of later calculations, we adopt s-wave separa- 
ble potential F(k— k') = C/ii;(k)u'(k'), and approxi- 
mate u;(k) = 1/Vl + e"(|k|-'=o)/fco with a > 1. The 
s-wave scattering length Og is related to U as m/(47ras) = 
1/U + {1/n) Efk|=o V(2ek), where fc, = fco ln(l + e")/a. 
For a repulsive interaction J7 > 0, Os is positive but up- 
per bounded by TT/(2kc) at C/ — > +00, and no bound 
state exists; for an attractive C/ < 0, diverges at Uc = 
—2Tr'^/{mkc), and only a sufficient attraction U < Uc 
with flg > can support a two-body bound state. 

(Ill) Continuum model with zero-range interaction 
potential in three dimension. H = ekcj^o-'^ko- + 
sEq.k.k' Cq+kt^q-kiCq-k'iCq+k't, and g is related to 
by m/(47ra,) = 1/g + (1/17) Y.^^^o l/(2ek). There is al- 
ways a two-body bound state when Og > 0. 

The single-band Hubbard model, as a simplified model 
for cold atoms in optical lattices (valid when interaction 
smaller than band gap) and many correlated materials, 
has been extensively studied before[3, l9l-[TTI|. The com- 
parison to previous known results justifies the validity of 
our approach and calibrates the correlation incorporated 
in our variational wavefrmction (w.f.). Then we apply our 
method to continuum models which are commonly-used 
for quantum gases and are of our primary interests. 

The variational w.f. we used for one down-spin system 
is similar to that used in the discussion of imbalanced 
Fermi gases [11], which is 

I*) = ^OC^^l + E '^kq4,+q_kt4tCqt + • • • 1^) 

y k>kF.q<kF / 

(1) 

where \N) represents a Fermi sea of iV-spin up particles 
with Fermi momentum kp. ". . ." in ^ represents terms 



2 



contains more than one particle-hole pairs of spin-up par- 
ticles. We compute the energy E (measured from energy 
of lA'^)) of l^*), and compare it with i?F {— -E|Ar+i)— E'iat)). 
Our main results are summarized as follows. 

First, for the single-band Hubbard model, we show 
that under certain conditions, iJintl'I') = and £ — 
{^\Ht\^) < E-p for most range of particle filling, except 
for nearby half-filling where the ground state is rigor- 
ously proved to be ferromagnetic by Nakaoka [1]. (See 
Fig. 1). This result agrees with previous studies by 
various other methods |9l-[ll|: We show that the real 
space representation of \^) corresponds to Gutzwiller's 
w.f. with optimized "backflow" type corrections, which 
implies \^) includes short-range correlations at least as 
the Gutzwiller projection. Similarly for the finite-range 
continuum model, we find £ < E-p for all range of C/ > 0. 
(See Fig. 2). This shows that generally in a purely repul- 
sive interaction model, a fully magnetized ferromagnetic 
state can not be ground state even for infinite U , which 
is in sharp contrast to Stoner's mean-field conclusion. 

Second, for the zero-range continuum model, apart 
from the polaron state with negative energy discussed 
before in Ref. [l2| - [l5j . we find the w.f. orthogonal to 
the polaron state has a positive £ for as > 0. At small 
feptts, £ as a function of fepOs follows the prediction of 
perturbation expansion very well. At large fepOg, £ sat- 
urates (due to unitarity) to 1.82i?F- Below k-pa^ = 2.35, 
we can find variational state with which £ < Ep , but fail 
to find such cases otherwise. (See Fig. 3). This result 
contradicts to the previous prediction kpa^ = 1.40 based 
on perturbation theory 0, hj. Thus the necessary con- 
dition for ferromagnetism is a large enough > 2.S5/kp, 
which requires a sufficiently attractive interaction poten- 
tial to support a bound state and cause resonant scatter- 
ing, and more importantly the system has to be prepared 
in the metastable scattering state that is orthogonal to 
molecule state. Nevertheless, our method can not prove 
this is a sufficient condition. Similar result is found for 
the finite-range model with U < and nearby scattering 
resonance. The details of our calculation are explained 
below. 

Single-band Hubbard model: Here we can first prove 
Theorem 1: Under two conditions that (i) 0kq = 4>k 
independent of q, and (ii) (po = — X]k>kF "^k, the w.f. 



00 c, 



qoJ. 



E 

k>kF,q<kF 



is an exact eigenstate of Hmt, and -ffint|V')i = 0. 

This theorem can be verified straightforwardly. The 
subscript i of j?/') means it contains one particle-hole pairs 
of up-spins, and the best variational energy for this state 
is denoted by £^^\ Condition (i) ensures that Hint acting 
onto will not generate two particle- hole term, and 
condition (ii) ensures zero interaction energy. For an in- 
tuitive understanding of this zero interaction energy, we 
can Fourier transfer it to real space. Using condition (ii). 




FIG. 1: 5E = £^^'> - E^iin units of t) as a function of particle 
density p = N/Ns in a 3D cubic (solid) and a 2D square 
(dashed line) lattice. 



one can show the w.f. ^ is equivalent to 



Here Vm = 1 ~ '^m-t'"^"^'t presents standard Gutzwiller 
projection operator. Ns is the number of lattice sites, 
'/'mn = X]k>kF 4'k^^^^"^~'^\ ^ud the sccoud term presents 
"backflow" type corrections. It becomes obvious that 
there will be no double occupancy and no interaction 
energy. Hence, we have established a momentum space 
w.f. representation of short-range correlation, which can 
be generalized to free space straightforwardly. 

Next we try to find the minimum total(=kinetic) en- 
ergy by introducing a Lagrange multiplier that can be 
proved as just the energy £^^^ (counted from the band 
bottom), T = {iplHtl^} - £'^'^'>{ip\ip). Minimization of T 
gives qo = and the self-consistent equation 



E 

k>kF 



£(1) 



Eq<kF(^kq-fW; 



- 1, 



(3) 



in which ekq = fq-k + Ck ^ Cq- Its solution gives f 

SE = £^^^ — Ep as a function of particle density p in 
a cubic and a square lattice is plotted in Fig. [1] We 
find for p < pc, 5E < means the ferromagnetic state 
is always unstable even for infinite U. In fact, for a cu- 
bic lattice, when p — > 0, £'^^^ goes like p while Ep goes 
like p^/^, hence 6E is negative and shows a much rapid 
decrease compared with the square lattice. While for 
p > Pc 5E > 0, this is consistent with Nagaoka theorem, 
which forces SE to be positive when p — > 1. pc obtained 
as 0.59 (0.76) for the square (cubic) lattice agree with 
previous studies by a finite-size real space evaluation [loj 
or by using Green functions [TT|. By contrast, our varia- 
tional w.f. (Eq. [5]) greatly simplifies the calculation, and 
also enable calculations in the thermodynamic limit sim- 
ply by employing density of state in a one-dimensional 
integral equation. Moreover, this w.f. can be systemat- 
ically improved by including multiple particle-hole con- 
tributions. 

Theorem 2: Consider the wave function \'ip)n that con- 
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tains up to n particle-hole pairs 

E 



1.0 



m=0 ^ ' {ki}5";{qj};' 



/,(™) t 



with pm = qo+X]"=i(Qi~l^i)' where {ki}™ denotes a set 
{ki, . . . , km} and {qj}™ denotes {qi, . . . , qm}- It satis- 
fies i?int|V')n = under the condition (i) '/'{■k-\" {qj}'" 
be expressed as B 



can 



(m) 

{ki}5";q„_i, 



(m) 

{ki}™:q„,q„_2 



{ki}™;qm,qm-i, 



C 



S(o) = 0, = 



(m-1) 



(m) 



(n) 

{ki}r;{qj}r 



(m) 

{ki}r;{qj}r 



• ,qi 
with 



{ki}5''-\k„;{qj}7 



= 0; and (ii) 
1-1 — for any 



{ki}r-^{qj}r-^ ' 

{k,}T-' and {qj}!"-!. 

Similar to Theorem 1, this theorem can be verified 
straightforwardly, and condition (i) ensures H-mt acting 
on \4')n will not generate n + 1 particle-hole terms, and 
condition (ii) ensures zero-energy. Obviously, |V')n-i is 
a special case of \ip)n, and hence < f (""I), where 
f '^"^ is the best variational energy for \ip)n- By including 
more particle-hole terms, one can always further lower £, 
and make pc more close to unity. On the other hand, one 
can also prove limp_j.o(f — — )■ 0, namely, 
the multiple particle-hole contribution vanishes and 
becomes exact at low density limit. 

Finite-range continuum model: We consider the same 
variational w.f. as Eq. [2] with 0kq = 0k- The saddle 
point equation from energy minimization is 

^Wk E ( '?^o + E ^k'0k' I = - ekq)(/>k, (4) 

q \ k' / q 

where we have adopted the approximation Wq ~ 1 for 
|q| < Ikp] < fco. From Eqs. 0] one can obtain a self- 
consistent equation 



2Ef 



1 



+ 00 



P-£/i2EF)' 



(5) 



where k = k/k-p and Uq — Um,k-p / {2Tr'^) are dimension- 
less. The solution to Eq. [5] is plotted in Fig. [2ja). Fig. 
^h) shows the results are insensitive to the choice of 
the free parameter a in the separable potential. As one 
can see, for Uq — > +oo^ £ saturates to a finite value, 
which is smaller than E^. This saturation precisely re- 
fiects the physics that the repulsive interaction can be 
strongly renormalized by short-range correlations, whose 
upper bound is the kinetic energy cost for the screen- 
ing, and is always finite even when the bare interaction 
diverges. In fact, for /eq 3> k-p, the integral at r.h.s. 
of Eq. [5] is approximately k^/k-p, hence at f7 — > +00, 
£ / Ep — 2kp/{3ko) < 1. For a general case, we show in 
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FIG. 2: (color online). £/Ey for the continuum model with 
finite-range interactions, (a) £^ as a function of Uo for ko/kp = 
2 or 5, a = 10. Arrows denote the corresponding critical Uc 
that can afford a two-body bound state, (b) the same plot as 
(a) for fco/fcp = 5 and q = 2 or 10. (c) £/Ep at Uo +00 as 
a functions fco for q = 2, 3, 10. 



Fig- IHc) £/Ep at Uq — +00 as a function of k^/kp at 
different a, and this ratio is always below unity. 

Zero-range continuum model: This model is in fact the 
/co +00 limit of the finite-range model with U < Q. 
For the w.f. of Eq. [21 = flV') gives the same self- 

consistent equation as in Ref. [l2|. 



_£_ 

'e^ 



dqq 



TT 1 



£ 



2 kpUs V^F 

where g{£ / Ep, q) is given by the integral 

P sin 9 



(6) 



dk / de^ 

Wo 2k^ -2qkcose-£/Ep 



where q — q/kp and k = k/kp. For this model there is 
always a bound state in a two-body problem for as > 0, 
and correspondingly, there is always a polaron solution 
with negative £ when as > 0, which has been exten- 
sively discussed in Ref. |12h15I |. In addition, there is 
always a positive energy solution orthogonal to the po- 
laron solution. For small kpcis applying the second- 
order perturbation theory 16] to this case, one will obtain 
£ /Ep = Akpas/ (Stt) -f 2{kpasY / ir"^ . It is found our varia- 
tional results fit quite well to this expansion in the regime 



\k 



pOs 



< 



0.6 (see inset of Fig. [3]), and start to deviate sub- 
stantially when kptts ^ 0.6. At fepOs — ^ -l-oo, we find the 
energy under this trial w.f. saturates at 1.82£'f instead 
of diverging as perturbation theory predicts. This sat- 
uration is due to unitary limit of resonance interaction. 
We find £ < Ep for kpUg < 2.35, which means a fully 
polarized ferromagnetism is definitely unstable below a 
critical kpa^ = 2.35. It contradicts to previous results 
based on a second-order perturbation ^7] which predicts 
the system becomes fully magnetized at fcpas = 1.40 [13]- 
The reason for this discrepancy is because the second or- 
der perturbation overestimates the interaction effects in 
the regime fcpas > 1- Similar behavior is found for a res- 
onance scattering with finite interaction range, as shown 
in the Uo < part of Fig. ^a) . 
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FIG. 3: £/Ef as a function of kpa^ for the continuum model 
with zero range interactions. Dashed horizontal lines denote 
saturated values of £/Ef at unitarity. Inset is fit to the per- 
turbation results (see text) for small fepOs. 

We would like to point out several intrinsic relations 
between (a) single-band Hubbard model at [/ ^ t, (b) 
finite-range continuum model at Uq = +oo and (c) zero- 
range continuum model at = -|-oo. First is the relation 
between the coefficient of different terms in the w. f. 
(Eq. [2]): for (b) we have (f>o + X]k>kF ^k^k = 0; this 
is equivalent to (a) with all Wk — 1 and a momentum 
cut-off imposed by tt/ol instead of ko (ol is the lattice 
constant); for (c) we can define — 0o ~^ ^!^k>kF 
then we have (pkq = Xq/(^ ~ ^kq) with Xq = gXq/^^, 
so at unitary limit one can get similar relation as (/>o + 
Ek>kp(<^kq + X°/(2ek)) + Eq<kp l/(2eq) = 0. Sccoud, 
for (a) at low density limit and (b) with kp ^ ko we 
can show £/Ep ~ kpTeS, here the effective range Toff — 
ol for (a), and l/fco for (b), which directly leads to the 
instability of ferromagnetism in these limits; while (c) 



does not fall into the class of pure repulsive interactions. 

The calculation above uses the w.f. that contains only 
one particle-hole pairs of up-spins. Including multiple 
particle-hole pairs can systematically improve the results, 
lead to further lower £ and increased fepOg. However, 
from the experience in studies of polaron branch, it is 
found the single particle-hole w.f. can already produce a 
result very close to Monte Carlo simulation [14| and later 
experiments p^ . and the reason for this perfect agree- 
ment is understood as a nearly perfect destructive inter- 
ference of higher-order particle-hole contributions [isj . 
Hence it is unlikely £/Ep at resonance can be reduced 
from 1.82 to below unity. Even though, we only prove 
ferromagnetic is stable against single spin flip above a 
critical kpa^, and we can not rule out the possibility that 
a state with more down-spin can be energetically more 
favorable due to more exotic correlations. 

Our results bring forward two intriguing issues, (i) for 
0-6 ^ kpGs < 2.35, the system can neither be fully ferro- 
magnetic nor be well described by perturbation theory. 
It is in a very interesting strongly interacting quantum 
phase with large ferromagnetic and/or short-range fluc- 
tuations; (ii) for kptts > 2.35, our approach predicts the 
system will become ferromagnetic. If it is true, what is 
the smoking gun experimental evidence? If future ex- 
periments flnd it is not true, then it means the system 
contains much stronger correlations than discussed here. 
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